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Abstract 

We first analyse the restricted four-body problem consisting of the Earth, the Moon and the 
Sun as the primaries and a spacecraft as the planetoid. This scheme allows us to take into account 
the solar perturbation in the description of the motion of a spacecraft in the vicinity of the stable 
Earth-Moon libration points L 4 and L 5 both in the classical regime and in the context of effective 
field theories of gravity. A vehicle initially placed at L 4 or L 5 will not remain near the respective 
points. In particular, in the classical case the vehicle moves on a trajectory about the libration 
points for at least 700 days before escaping away. We show that this is true also if the modified 
long-distance Newtonian potential of effective gravity is employed. We also evaluate the impulse 
required to cancel out the perturbing force due to the Sun in order to force the spacecraft to 
stay precisely at L 4 or L 5 . It turns out that this value is slightly modified with respect to the 
corresponding Newtonian one. In the second part of the paper, we first evaluate the location of 
all Lagrangian points in the Earth-Moon system within the framework of general relativity. For 
the points L 4 and L 5 , the corrections of coordinates are of order a few millimeters and describe 
a tiny departure from the equilateral triangle. After that, we set up a scheme where the theory 
which is quantum corrected has as its classical counterpart the Einstein theory, instead of the 
Newtonian one. In other words, we deal with a theory involving quantum corrections to Einstein 
gravity, rather than to Newtonian gravity. By virtue of the effective-gravity correction to the long¬ 
distance form of the potential among two point masses, all terms involving the ratio between the 
gravitational radius of the primary and its separation from the planetoid get modified. Within 
this framework, for the Lagrangian points of stable equilibrium, we find quantum corrections of 
order two millimeters, whereas for Lagrangian points of unstable equilibrium we find quantum 
corrections below a millimeter. Finally, general relativity corrections to Newtonian position of 
collinear Lagrangian points turn out to be below the millimiter, whereas on stable equilibrium 
points they are of order of a few millimiters. 

PACS numbers: 04.60.Ds, 95.10.Ce 
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I. INTRODUCTION 


In the space surrounding two bodies that orbit about their mutual mass center there are 
five points where a third body will remain in equilibrium under the gravitational attraction 
of the other two bodies. These points are called Lagrangian points in honour of Joseph 
Lagrange, who discovered them in 1772 while studying the restricted problem formed by 
the Sun-Jupiter system. The discovery of their physical realization, i.e. the Trojan group 
of asteroids, began only in 1906 thanks to the astronomer Max Wolf with the first-seen 
member of this group, 588 Achilles, which is located near the triangular libration point 
of the Sun-Jupiter system. Today we know that there are 3898 known Trojans at the 
triangular Lagrangian point L 4 and 2049 at L 5 [1]. In the sixties, simultaneously with the 
increased interest in space explorations, the question of existence of Lagrangian points with 
respect to other primaries, especially for the Earth-Moon system, arose quite naturally. In 
fact, if there are stable stationary solutions for various primary combinations, then from 
a practical point of view placing observational platforms at these points becomes feasible, 
especially in a really close and accessible system like the Earth-Moon system, which is also 
the most convenient system from an economic point of view. While the Sun-Jupiter system 
clearly possesses a collection of asteroids at the triangular libration points, the ability of 
the Earth-Moon system to collect debris or dust at the corresponding points and in what 
is called Kordylewski clouds is still in question (see Ref. j^J] for further details). The major 
perturbing effect on the Trojans is represented by Saturn, while the stabilizing forces come 
from the Sun and Jupiter. The major perturbation on the Earth-Moon libration clouds is 
the Sun and the stabilizing effects are derived from the Earth and the Moon. This explains 
why the existence of accumulated material at L 4 or L 5 in the Earth-Moon system is not 
so obvious. Bodies at the triangular libration points of the system consisting of the Sun 
and another planet would face the perturbations from Jupiter; therefore, it is not surprising 
that the only currently known material accumulation is confined to the Sun-Jupiter system, 
although some asteroids were found also in the Sun-Earth system around the libration point 
L 4 , as is shown by recent observations Q. As far as the collincar Lagrangian points for the 
Earth-Moon system are concerned, we know that L\ allows comparatively easy access to 
Lunar and Earth orbits with minimal change in velocity and has this as an advantage to 
position a half-way manned space station intended to help transport cargo and personnel 
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to the Moon and backwards, whereas L 2 would be a good location for a communications 
satellite covering the Moon’s far side and would be an ideal location for a propellant depot 


as part of the proposed depot-based space transportation architecture 




a 


of US 


Recently, inspired by the works in Refs. |5|-|l2] on effective held theories of gravity, some 
ied this theoretical analysis to the macroscopic bodies occurring in 
is| . especially in the Earth-Moon system. It has been demonstrated 
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15| have app 


celestial mechanics 
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that in the quantum regime, when only the interaction potential is modified in the La- 
grangian of Newtonian gravity, the position of collinear Lagrangian points is governed by 
four algebraic ninth degree equations, which reduce to two algebraic fifth degree equations 
in the classical regime, while the quantum corrected position of the noncollinear libration 
points is described in terms of a pair of quintic equations, which predict that the classical 
equilateral triangle picture is no longer valid in the quantum scheme. For the Earth-Moon 
system, the prediction about the discrepancy between classical and quantum corrected quan¬ 
tities is of the order of millimeters. This magnitude is comparable with the instrumental 
accuracy of point-to-point laser Time-of-Flight (ToF) measurements in space typical of the 
modern Satellite/Lunar Laser Ranging (SLR/LLR) techniques [19h32]). The full positioning 
error budget of the orbits of satellites equipped with laser retroreflectors and reconstructed 
by laser ranging depends also on other sources of uncertainty (related to the specific orbit, 
satellite and retroreflector array), in addition to the pure point-to-point laser ToF instru¬ 
mental accuracy (related to the network of laser ranging ground stations of the ILRS, i.e., 
the International Laser Ranging Service j^oj] . The full positioning error budget can be larger 
than millimeters. 

This is an interesting potentiality, because we are dealing with predictions which might 
become testable in the Earth-Moon system. This is a novel feature in the theory of quantum 


gravity, because all other theories are so far unable to produce testable effects |33l-l43j. These 


predictions become more realistic if we include the perturbations due to the gravitational 
presence of the Sun, in other words we have to face up the restricted problem of four bodies, 
consisting of the Sun, the Earth and the Moon as the three primaries and the fourth body 


e.g. a laser-ranged test mass, a spacecraft or by exploiting the solar sail technology 


15 


44 


5l|]) which has an infinitesimal mass, to avoid affecting the motion of the primaries. As 


we know, in the restricted three-body problem the motion of the two primaries is exactly 
described by the equations of motion governing the two-body problem. Therefore, we may 
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generalize the problem first by solving the dynamical equations describing the motion of the 
three primaries and then by finding the motion of the planetoid in the presumably known 
gravitational held produced by the primaries. Since no closed-form solution is known for the 
full three-body problem, this generalization to the case of four bodies is rather difficult. A 
practicable possibility consists in assuming the motions of the three primaries and, without 
attempting to establish the exact solution of the equations governing these motions, accept 
an approximate solution. Such an approximation may be, for instance, that the Earth and 
the Moon move in elliptic orbits around their mass center and that the mass center of the 
Earth-Moon system, in turn, moves in elliptic orbit around the Sun. The plane of the 
orbit of the mass center of the Earth-Moon system, which is called the plane of ecliptic, is 
inclined relative to the plane containing the orbits of the Earth and the Moon. A simpler 
approximation would consist in neglecting the eccentricity of all orbits, i.e. assuming that the 
Earth, the Moon and their mass center have circular orbits. Under these assumptions, the 
authors of Ref. [52] did show that, although it is widely accepted that, with the introduction 
of the Sun, the points L 4 and L 5 of the Earth-Moon system cease to be equilibrium points, 
stable motion may be possible in a region around these noncollinear libration points. The 
term “stable” here indicates that the planetoid will remain within a certain region for the 
period of time during which the motion is studied. The work in Ref. [52J demonstrated that 
a spacecraft moves on a trajectory around the stable libration points for at least 700 days 
before the solar influence causes it to move through wide departures from the Lagrangian 
points. Indeed, from the analysis of the plots it does not appear that, after 700 days, a 
limiting value for the envelope is approached. It would be interesting to recover this feature 
directly from the solution of the dynamical equations (if they were known), since at the 
present state we believe that the form of the equations involved (see Sec. II) does not allow, 
by itself, such a deduction. 

The first purpose of our paper consists in showing that this is true also if we assume the 
quantum corrected potential discovered in Refs. 5|, LL 2 ||, and Secs. II and III are devoted to 
this topic. All the considerations made in these Sections represent the natural extension of 
our previous papers, as we continue to describe the three-body problem in the context of 
effective held theories of gravity by adding all features that would contribute to make this 
subject as close as possible to reality, in order to encourage the launch of future space mis¬ 
sions that could verify the model we are proposing. On the other hand, one has to consider 
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that general relativity is currently the most successful gravitational theory describing the 
nature of space and time, and well confirmed by observations. In fact, it has been brightly 
confirmed by all the so-called “classical” tests, i.e. the perihelion shift of Mercury, the de¬ 
flection of light and the Shapiro time delay, and it has also gone through the systematic test 
offered by the binary pulsar system “PSR 1913+16”, since the orbit decay of this system is 
perfectly in accordance with the theoretical decay due to the emission of gravitational waves, 


as predicted by general relativity. Furthermore^ 


renewed interests for relativistic astrophysics 


53 


agrangian points have recently attracted 


56], where the position and the stability 


of Lagrangian points is described within the post-Newtonian regime. For all these reasons, 
we believe that our model is incomplete without a comparison with the Einstein theory. 
Therefore, Sec. IV studies all Lagrangian points within the framework of general relativity, 
to establish the most accurate classical counterpart of the putative quantum framework that 
we have set up. By taking seriously into account the important role played by the Einstein 
theory within this scheme, in the last part of this paper we describe a new quantum corrected 


regime where the underlying classical theory is representec 
Newtonian theory. All the considerations made in Refs. 


by 


general relativity, rather than 


13Ml5l]. in fact, are characterized 


by the fact that the classical theory for which quantum corrections are computed is the 
Newtonian theory, instead of Einstein’s one. But, if general relativity is the most successful 
classical theory of gravitation, then we have to consider a scheme where quantum correc¬ 
tions to general relativity are evaluated. This topic is investigated in Sec. V, where we also 
show that, among all quantum coefficients K]and Ko in the long-distance corrections to the 
Newtonian potential available in literature |5|, 1214151]. the most suitable ones to describe 
the gravitational interactions involving (at least) three bodies are those connected to the 
bound-states potential. Finally, conclusions and open problems are discussed in Sec. VI. 


II. THE QUANTUM CORRECTED EQUATIONS OF THE RESTRICTED FOUR- 
BODY PROBLEM 


We start by introducing the classical dynamical equations governing: the motion of the 
planetoid in the gravitational held of the Earth, the Moon and the Sun [52|, |57[ • We suppose 
that the Earth and the Moon move in circular orbit around their mass center and the mass 
center, in turn, moves in circular orbit about the Sun. The Earth-Moon orbit plane is inclined 
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at an angle i = 5°9' to the plane of the ecliptic. We introduce the rotating coordinate system 
£, 77 , £ with the Earth-Moon mass center as its origin and characterized by the fact that the 
£ axis lies along the Earth-Moon line, the 77 axis lies in the Earth-Moon orbit plane and the 
£ axis points in the direction of the angular velocity vector of the Earth-Moon configuration. 
The £, 77 axes rotate about the £ axis with the angular velocity u: of the Earth-Moon line. 
If the vector 1Z = (£, 77 , £) indicates in this coordinate system the position of a spacecraft of 
infinitesimal mass, the vector dynamical equation describing its motion is 

77-1 - oj x (2 7?. + w x 77) = —V^E -I- \ZtzU 4 - S, (2-1) 


where 


Gmi Gni2 
Pi P 2 


U = G 1713 


1 

Ps 


7Z ■ TZj, 


( 2 . 2 ) 

(2.3) 


with G being the universal gravitation constant, mi,m 2 and m 3 the mass of the Earth, the 
Moon and the Sun, respectively, pi, p 2 and p 3 the distances from the planetoid of the Earth, 
the Moon and the Sun, respectively, 77 3 the distance of the Sun from the Earth-Moon mass 
center, and lastly S describes the solar radiation pressure. Written in components, Eq. (12.ip 
becomes 


£ - 2 ufj - w 2 £ 
77 + 2cn£ — u 2 ri 

C 


dv du „ 

-~~ds + a? + Ss ’ 

_dV_ dU_ 
dp dr/ Vl 
dV dU n 

— 7TT + "TTV + *Sc- 

d( d( 


(2.4) 

(2.5) 

( 2 . 6 ) 


We can write Eqs (I2.4p ~ 02.6p in what we denote by x, y, z system, which is the rotating 
noninertial coordinate frame of reference centered at one of the two noncollinear Lagrangian 
points, e.g. L 4 . If we use the transformations 


s T £pi 


v = y + Vp , 
C = z, 
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(2.7) 













where £ p and y p are the constant coordinates of the libration point L 4 in the £, rj, ( system, 
then Eqs. (12.4ft - (12.61) become 


Gm- 

x = 2 uy + (x + £ p )oj 2 - (x 3 + £ p )(hL) 2 + S x + ^2 — _ x )i (2.8) 

i =1 7k 

3 ^ Qjji- 

y = -2ux + (y + p p )uj 2 - ( 2/3 + p p )(f0 ,) 2 + S v + —^(j/i - I/), (2.9) 

»=1 Pi 

3 ^ 

2 = -* 3 (SJ„) 2 + & + V -P(z, - 2 ), (2.10) 

tr « 

where is the angular velocity of the Earth-Moon mass center around the Sun, and the 
relation Gm 3 / (JZ 3 ) ? ' = (hL ) 2 has been exploited. Moreover, the distances p,. are given by 

(pi) 2 = (xi - x) 2 + fa - y) 2 + (zi - zf (i = 1, 2, 3), (2.11) 

where the coordinates (xi,yi) and (a: 2 , 2 / 2 ) of the Earth and the Moon respectively are 
deduced from (12. 7 j) once the coordinates (£ p , r] p ) of L 4 are known (remember we have z\ = 
z -2 — 0), whereas the coordinates of the Sun are given by the relations 

x 3 = 7^3 ( cos 0 cos $ + cos i sin 0 sin 6 ) — £ p , 

2/3 = — Tl 3 (cos 0 sin 9 — cos i sin 0 cos 9) — rj p , (2-12) 

z 3 — 703 sin 0 sin i, 


where 0 is the angular position of the Sun with respect to the vernal equinox and measured 
in the plane of the ecliptic, and 6 describes the position of the Earth-Moon line with respect 
to the vernal equinox measured in the Earth-Moon orbit plane (see Fig. 2 of Ref. 52]). The 
relations defining these angles are 


0 = f00 + 00, 

0 = + 6 q, 


(2.13) 


where 0o and do are the initial values of 6 and 0. For our computation we have used the 
following numerical values: = 1.99082 x 10 ~ 7 rad/s, ui = 2.665075637 x 10 -6 rad/s, 

0o = p o — 0 (i.e. the initial position of the Sun will be on the extended Earth-Moon line, 
with the Moon in between Earth and Sun). Moreover, following Ref. 15|, we have the 
classical values £ p = 1.87528148802 x 10 8 m and p p = 3.32900165215 x 10 8 m. If we initially 
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set S' = 0 in dm we obtain that the perturbative effect of the Sun makes the spacecraft 
ultimately escape from the stable equilibrium point after about 700 days j^], as is shown 
in Figs, m and 0 As can be noticed from Fig. [Q the irregular initial motion damps out 



x 


FIG. 1: Parametric plot of the spacecraft motion about L 4 resulting from zero initial 
displacement and velocity in the classical case. The quantities appearing on the axes are 
measured in meters and the time interval considered is about 4 x 10 7 s. 



FIG. 2: Plot of the spacecraft motion about L 4 in the ^-direction resulting from zero initial 
displacement and velocity in the classical case. The quantities on the axes are measured in 

meters and in seconds. 

and there is an approximate one-month periodicity associated with the motion. Moreover, 
Fig. [2] shows that the amplitude of the motion increases with time and that the period of 
motion is about 27.6 days, a value really near to the 29.53 days of the synodical month. 
All these results indicate that the spacecraft will escape from the equilibrium point L 4 (or 
equivalently L 5 ) or, in other words, the perturbing presence of the Sun makes the points L 4 
and L 5 cease to be equilibrium points, but they are “stable” in the sense indicated in the 
Introduction. 

All these considerations are valid within the classical scheme, whereas in the quantum 
corrected regime (see Secs. IV and V) the Newtonian potential is corrected by a Poincare 
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asymptotic expansion involving integer powers of G only, so that Eq. ( 12 .ip can be replaced 
by the vector dynamical equation 


77-1- cu x (2 77 -I - oj x 7?.) — —-I- V pJJq + S , 


(2.14) 


with 
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1 ^1 
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k 2 

1 Gm 2 

r k\ k 2 

1 + + 

Pi 

L pi 

(piA 

P2 

L P 2 (p 2 ) 2 J 


U q = 


Gm 3 

P3 


k'[ 

1 + — + 
P3 


(Ps) 2 , 


Gm. 


77 • 7^3 

'W 


2k'; 

1 + + 


3 k 2 


K 3 ' (77 3 ) 2 , 


(2.15) 

(2.16) 


where, following Ref. [12j|, we decide to adopt the results concerning the bound-states 


potential. Even without knowing the detailed calculations of Sec. V, we may point out 


that, in classical gravity, the Levi-Civita cancellation theorem 


holds, according to which 


the IV-body Lagrangian in general relativity can be always reduced to a Lagrangian of N 
material points. In other words, it is not necessary to assume that we deal with point 
particles for simplicity, but the effects of their size get eventually and exactly cancelled. 
Now the quantum corrections considered in Refs. 


13 


15| deal precisely with the long¬ 


distance Newtonian potential among two such masses, and consider three distinct physical 
settings: scattering, or bound states, or one-particle reducible [l2[. We think that, in celestial 
mechanics, the bound states picture is more appropriate for studying stable and unstable 
equilibrium points. Therefore we set (cf. Sec. V) 


h = 


Gmi . Gm 2 


- k' — 

1 "a 


, — 


Gm 3 41 2 

k‘2 = —(Ip) , 


(2.17) 


2 c 2 ’ 1 2 c 2 ’ 1 2 c 2 * ■'" IOtt ' 

lp being the Planck length. The occurrence of the term fc 2 , which is quadratic in the 
Planck length, cannot be obvious for the general reader, and hence we here summarize its 
properties and derivation, following our sources js-7, 12]. The one-loop quantum correction 
to the gravitational potential is a low-energy property independent of the ultimate high- 
energy theory. The potential of gravitational scattering of two heavy masses turns out to 
be 


V(r) = 


GMm 


1 + 3 


G(M 


m. 


c 2 r 


+ kl 

o 


From dimensional analysis one can indeed expect a term like %, because the unique dimen¬ 
sionless term linear in K and linear in G is The classical post-Newtonian correction is 
also a well-known dimensionless combination, without h. We have a numerical factor of —5 
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in (2.17) obtained as 3 — ^ = — A, because the bound-state contribution 12 ] written within 
square brackets above is — By Fourier transform, the corresponding results in momentum 
space turn out to be [7] 

- ->• 4> 4^4 x 4^4 x ? 2io g(? 2 )’ ^ 3 (^) ^ 4 x ^ 2 - 

r q r q 2 r s q q 2 

The one-loop potential is obtained from one-graviton exchange, with the dy resulting from 
the massless propagator. The corrections linear in h are due to all one-loop diagrams that 
can contribute to the scattering of two masses. The kinematic dependence of the loops then 
brings in nonanalytic corrections of the form Gm^/q 2 , Gq 2 log(g 2 ), as well as analytic terms 
Gq 2 . ffowever, the Fourier transform of the analytic term is a Dirac delta in position space, 
and hence analytic terms do not contribute to long-distance modifications of the potential. 


The above correspondences are made precise by the following integrals 


12 ]: 


d 3 g 


qqr __ 


d 3 g 


q q-r __ 


d 3 g 


e^logd q\ 2 ) = 


J (27r) 3 \q\ 2 47rr’ J (27r) 3 |g| 27r 2 r 2 ’ J (27r) 3 27Tr 3 

In the x,y,z system, instead of Eqs. fl2.8h - fl2.10h . Eq. (12.1 4j) . written in components, 
gives rise to the system 


x = 2uy F (x F ^p)cj 2 - (x 3 F £ p )(hL) s 
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P3 (Ps) 2 , 


+ S'z, 


( 2 . 20 ) 


where we have used the fact that z\ = z 2 = 0. Setting S' = 0, we have integrated Eqs. (12.18h — 
(I2.20h and we have discovered that the situation is almost the same as in the classical case 
(see Figs. [3] and |4]) , i.e. the planetoid is destined to run away from the triangular libration 
points in about 700 days. This means that, also within a quantum corrected scheme, the 
gravitational effect of the Sun spoils the equilibrium condition at L 4 and L 5 . 
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FIG. 3: Parametric plot of the spacecraft motion about L 4 resulting from zero initial 
displacement and velocity in the quantum case. The quantities appearing on the axes are 
measured in meters and the time interval considered is about 4 x 10 7 s. 



FIG. 4: Plot of the spacecraft motion about L 4 in the ^-direction resulting from zero initial 
displacement and velocity in the quantum case. The quantities on the axes are measured 

in meters and in seconds. 

III. THE SOLAR RADIATION PRESSURE AND THE LINEAR STABILITY AT 

L 4 


At this stage, we assume the presence of the radiation pressure both in the classical 
equations fl2.8p - fl2.10p and in the quantum ones fl2.18IMI2.20lh The solar radiation pressure 
is given by 


A 


s = -K—— jp 3 , (3.1) 

mypzY 

where A is the cross-sectional area normal to p^, m is the planetoid mass and K is a constant. 
Inspired by Ref. 52], we use the value K = 2,048936 x 10 17 1V. We have integrated the 
classical equations fl2.8p - fl2.10p and we have found that the presence of the solar radiation 
pressure causes the vehicle to move further away from L 4 in a given time, as one can see 
from Fig. [5] In particular, the larger the ratio A/m is, the larger the envelope of the motion 
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FIG. 5: Parametric plot in the classical regime of the spacecraft motion about L 4 in the 
presence of the solar radiation pressure and considering A/m = 0,159 m 2 /Kg. The initial 
displacement and velocity are zero. The quantities appearing on the axes are measured in 
meters and the time interval considered is about 1 x 10 7 s. 



FIG. 6 : Parametric plot in the quantum regime of the spacecraft motion about L 4 in the 
presence of the solar radiation pressure and considering A/m = 0,159 m 2 /Kg. The initial 
displacement and velocity are zero. The quantities appearing on the axes are measured in 
meters and the time interval considered is about 1 x 10 7 s. 


is 


52]. 


Interestingly, in the quantum case ruled by effective gravity the situation is a little bit 
different. Unlike the classical regime, the presence of the solar radiation pressure in the 
equations (I2.18jf (I2.20j) does not show itself through the fact that the spacecraft goes away 
from the triangular libration points more rapidly, but it results in a less chaotic and irregular 
motion about L 4 , which ultimately make the planetoid escape from L 4 , like in the classical 
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FIG. 7: Fig. 7a: Spacecraft motion about L 4 in the classical regime and with an initial 
velocity of 3 m/s at 60°; Fig. 7b: Spacecraft motion about L 4 in the classical regime and 
with an initial velocity of 3 m/s at 150°; Fig. 7c: Spacecraft motion about L 4 in the 
classical regime and with an initial velocity of 3 m/s at 240°; Fig. 7d: Spacecraft motion 
about L 4 in the classical regime and with an initial velocity of 3 m/s at 330°. 


case. These effects are clearly visible from Fig. Ef. We can also try to find the best set of 
initial conditions which leads to the smallest envelope of the motion of the planetoid. We 


have studied several sets of initial conditions both in the classical case anc 
one. In the classical regime, we completely agree with the results of Ref. 


in the quantum 


52j . We have found 


in fact that the amplitude of the spacecraft’s motion depends strongly on the position of the 
Sun (i.e. on the values assumed by 9o and i/jo) and on its initial position and velocity. For 
example, Fig. [7] shows the motion resulting from an initial zero displacement and different 


1 The different scale adopted in Fig. [G]with respect to the one of Fig. fallows us to better appreciate its 
features. 
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initial velocity ( 6 q = ^0 = 0) and the time interval considered is about 1 x 10 7 s. As we can 
see, the envelope of the motion in Fig. 7b is smaller at any time than the envelope of the 
motion shown in Fig. E3 

The situation is fairly the same in the quantum regime (Fig. [SJ), where we have discovered 
that one set of initial conditions (having 0 O = V’o — 0 ) exists which results in a smaller 
envelope of the spacecraft motion at any given time, as one can see from Fig. 8 b. This fact 
can be understood with a comparison between Figs. 6 and 8 b. The interesting difference 
with respect to the classical case consists in the fact that the reduction of the envelope of 
the planetoid motion produced by a nonzero initial velocity becomes more evident in the 
quantum regime. By inspection of Figs. 7 and 8 we discover a strong dependence on the 
initial conditions of the planetoid trajectories both in the classical and quantum regime. 
This suggests that, from an experimental point of view, it might be useful to drop off two 
or more satellites close to the Lagrangian points L 4 and L 5 with slightly different initial 
conditions for position and velocity. Measurements of the satellite differential positions, 
together with the measurement of the single orbits, could make it possible to discriminate 
between classical and quantum regime, without depending on the absolute knowledge of 
Lagrangian points’ location. 

If we want to force the particle to stay precisely at L 4 , we have to set aside the perturbing 
force due to the Sun by the application of a continuous force (see Fig. ED. Therefore, we 
have to study the following stability equation (in the £, 77 , ( system): 

- V n V + WtzU + S+- = 0, (3.2) 

m 

which becomes in the quantum case 

- V n V q + Vt z U q + A + — = 0, (3.3) 

m 

where m is the mass of the planetoid and F (respectively, F q ) represents the force to be 
applied to the spacecraft in order to make it stay precisely at L 4 in the classical (respectively, 
quantum) regime. If we consider Eqs. (13.2^ — fj3.3jl in the x, y, z coordinate system, we can 
exploit the simplification resulting from the fact that the planetoid must be at the position 
x = y = z = 0, hence Eq. (13. 2 j) . written in components, becomes 

£ p u ; 2 — (X 3 + £ p )(h2 ^) 2 + S x + H—- = 0, (3.4) 

^ J P; m 

i=1 
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FIG. 8: Fig. 8a: Spacecraft motion about L 4 in the quantum regime and with an initial 
velocity of 3 m/s at 60°; Fig. 8b: Spacecraft motion about L 4 in the quantum regime and 
with an initial velocity of 3 m/s at 150°; Fig. 8c: Spacecraft motion about L 4 in the 
quantum regime and with an initial velocity of 3 m/s at 240°; Fig. 8d: Spacecraft motion 
about L 4 in the quantum regime and with an initial velocity of 3 m/s at 330°. 




FIG. 9: Fig. 9a: Force per unit mass as a function of time required to induce stability at 
L 4 in the classical regime; Fig. 9b: Force per unit mass as a function of time required to 

induce stability at L 4 in the quantum regime. 
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13 C Wr) . 7? 

VpU 2 _ (?/3 + Vp)(^) 2 + Sy + - 3 —yi + — = 0, 

»=i 


m 


^(^oj) 2 + <Sz H-3—^3 H-— 0, 

p? m 


whereas from Eq. (13.31) we obtain 

C,pUJ (x 3 + £p)(f\;) 


2fc7 3k 2 
1 + -^ + 


Gm 2 x 2 

7FF 


l+ 2^ + 3 /c 2 


P2 


(P2) 2 


+ 


^3 ' (^ 3 ) 2 

Gm 3 x 3 
(p 3 ) 3 


+ 


Gmix 1 

“FF 

1+ 2^ + 3fc 2 


2£q 3fc 2 
1 + —+ 

Pi IPi 


(pO 2 


P3 


(P3) 2 , 


m 


(3.5) 

(3.6) 


+ S x + = 0, (3.7) 


r\ p u? - (y 3 + r] p )(fl u} ) 2 


2k" 

1 + — 
n 3 


3 k 2 

(FF. 


Gmiyi 

(Pi) 3 


2fc! 

1 + —+ 
Pi 


3fc 2 

FF, 


Gm 2 y 2 

(P 2) 3 


[1 + 1 

Gm 3 y 3 

2k" 

1H-- + 

3 k 2 

L P 2 (P 2 ) J 

(Ps) 3 

L Ps 

(Ps) 2 J 


F a 

+ Sy +^= 0, 

m 


- ^(^F 

Equations (13.41) 


2 k" 

1 + — 
^3 


3/c 2 

(F) 


+ 


2 k" 

1 + —+ 
P3 


3fc 2 

FF. 


+ ^ + ^ = 0. 
m 


(3.8) 


(3.9) 


and (13.71) 


Gm 3 z 3 
(P 3) 3 

make it possible for us to evaluate both the classical 


and the quantum force needed for stability and therefore the impulse per unit mass which 
the planetoid must be subjected to in order to stay in equilibrium exactly at L 4 . Bearing in 
mind that the impulse is defined as the integral of a force over the time interval for which 
it acts, and on considering a time interval of one year, we have found the following results 
for the classical and the quantum regime, respectively: 


I c /m = 747, 608255 N s/Kg, 
I q /m = 747, 608245 N s/Kg. 


(3.10) 


Of course, these considerations are preliminary because, even just at classical level, the four- 
body problem has only been studied from us within Newtonian gravity. We also note that 
this calculation suggests a gedanken experiment in which two satellites are sent to L 4 and L 5 , 
respectively. If the first satellite receives the impulse I c while the second receives the impulse 
I q , one might try to check, by direct comparison, which value is better suited for stabilizing 
the Lagrangian point, gaining support for classical or, instead, quantum theory. However, 
this configuration is merely ideal because, in light of the very small relative difference of 
the impulse in the two cases, it looks practically impossible to keep all the experimental 
conditions (satellite mass, actuator and readout calibration, initial conditions, solar radiation 
pressure etc.) identical within the required accuracy (less than 0.1 parts per million). 
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IV. THEORETICAL PREDICTIONS OF GENERAL RELATIVITY 


A. Noncollinear Lagrangian points 


The analysis of the previous section relies on the simple but nontrivial assumption that, 
since effective gravity modifies the long-distance Newtonian potential among bodies of 
masses and mg according to 2 

GmAmB 


V E (r) 

V E (r) 

c 2 m B 


r 


, , | (Ra + Rb) 

1 + I K\ -h K2 


o ipr 


+ 0(G 


Ra 

r 


1 + ( Kl 


(Ra + R B ) , ilp) 2 , 


T 


+ 0(G 


(4.1) 


for all values of r greater than a suitably large tq, where gravitational radii Ra, Rb and 
Planck length Ip are defined by 


_ Gm A D _ Gm B , _ 

ft A = ftB = 9—, tP = 


Gh 


(4.2) 


the resulting modification of Newtonian dynamics can be obtained by considering a classical 
Lagrangian where the Newtonian potential 

GmABlB 


V N (r ) = 


(4.3) 


is replaced by V E (r), while all other terms remain unaffected (cf. Sec. V). Although it would 
be inappropriate to use the quantum effective action to study the low-energy effects resulting 
from the asymptotic expansion (14.ip , the above assumption is a shortcut to describe a theory 
lying in between classical gravity and full quantum gravity. For this reason, it becomes 
important to study the predictions of classical gravity when general relativity is instead 


assumed. The work in Ref. 


551 ] has indeed done so by relying upon the Einstein-Infeld- 


Hoffmann equations of motion for a three-body system, but without studying an effective 
potential and the zeros of its gradient. However, such a potential is by now available in the 
literature, and the resulting approximate evaluation of Lagrangian points L 4 and L 5 was 
performed in Ref. [59|, while their stability in a suitable mass range was proved in Ref. 60]. 
Strictly speaking, in general relativity the libration points become quasi-libration points, 


2 We can say that the real number n\ is the effective-gravity weight of the sum of gravitational radii, whereas 
the real number «2 is the effective-gravity weight of Planck’s length squared. 
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61| for this feature, which reflects the expected emission of 


and we refer the reader to Ref. 
gravitational radiation. 

For our purposes, it is enough to consider the relativistic version of the circular restricted 
three-body problem in a plane, where, for primaries of masses a and (3 separated by a 
distance l, with gravitational radii R a = ?§- and Rp = and mass ratio p = | < 1, the 
classical angular frequency (or pulsation) u = is replaced by 59] 


n = u 


3 ( R a + Rp) 

~ 2 1 


i-i. 


P 


(4.4) 


3(1 + P) 2 

while, in a noninertial frame with origin at the mass center of the Earth-Moon system, the 
equations of motion of the planetoid in the planar case (with coordinates (£, rj)) read as 


£ — 2Vtr] = 
r) + 2h2£ = 


dW 

w 

dW 

dp 


d_ 
d t 

d 

d£ 


dW 

~d£ 

dW 

dfj 


(4.5) 

(4.6) 


where, upon denoting by r the distance of the planetoid from the primary of mass a (i.e. 
the Earth), and by s the distance of the planetoid from the primary of mass f3 (i.e. the 
Moon), given by 

2 


s = u 


e+ (l+J + "’ 

(4.7) 

-Ai )'♦* 

(4.8) 


one has the effective potential reading as [59] 


w = ^(( 2 + v 2 )+ C 2 


-Rq ^ -R/3 
r s 


[Ra? AM 

O ' O 




Rf) m(Arj + -n^ 


— + — ) f (€, v ,€, v ) 

r s 


(1 + P) 

Rp 


2/2 




r/ 


(1 + p) L 2 (1 + P) \ r 


P 


l 


(p-2) l + (l-2p)l 


rs 2(1 + p) r 2(1 + p) s 


where 


59] 


(4.9) 


(4.10) 


/(f.o.to) = r + >i 2 + 2Si(fi? - of) + fi 2 (f 2 + o 2 )- 

At all equilibrium points, the first and second time derivatives of coordinates (£, 77 ) should 
vanish, which implies that it is enough to evaluate the zeros of the gradient of VF(£, 77 ), 
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because 


59] 


d (dW 


d fdW 


= 0if£ = 7) = £ = p = 0. 


(4.11) 


dt \ d£ J d t \ dr) 

Note now that, by virtue of (14.7ft and (14.8ft . one has the formulas (IAlft - (lA4ft in the Appendix, 
and hence the two components of the gradient can be expressed in the form 


= Wi(f, 77 ,r) + W 2 (£,r),s) + Rpl ^ 


dW 

dr) 


= V 


(1 + p) rs [ \r 
W 3 (Z,r),r) + W 4 (f, 77 ,s) + 


n — — 

s ( _2 e 2 


l 


p 


RJ 3 q 2 


(1 + p) \r 2 s 2 


1 1 
+ 


(1 + p) rs \r 2 s 2 


, (4.12) 


(4.13) 


where the functions Wi ,..., W 4 are defined in Eqs. 


13 


. Thus, unlike the case of Refs. 
15], when the gradient of w is set to zero with t) ^ 0, one does not get an algebraic 


equation for r only. Since we are interested in numerical solutions of such an enlarged 
algebraic system with (at least) ten decimal digits, we set r = 7 /, s = FI, and we study the 
coupled algebraic equations for the real numbers 7 and T obtained from 


ynf = 0 , 

= 0 . 

r) or) 


(4.14) 


(4.15) 


where the fifth powers of 7 and T are suggested by the occurrence of terms proportional 
to 7 -5 and T " 5 in the derivatives and We can write Eqs. (I4.14p ~ 04.15p in a 


dr/ 


concise way, i.e. 


7 5 r 5 
7 5 r 5 


dW 


d£ 

1 dW 


EA.(rV = 0, j €{0,1,2, 3,4, 5}, 


n =0 
5 


77 dr) 


J2 B n(r j h n = 0, je{ 0,1, 2, 3,4, 5}, 


more 

(4.16) 

(4.17) 


n =0 


where the coefficients A n (T->) are given by 

fl 2 


A 5 (p) = r 5 


1 + 


2 c 2 


W + O 


£fi z + r 


3^ 7 

1 2 (1 + p) 


RpFl 2 


if ! 2 

2 


r 


cR 


l 2 


7 li l 2 ( 2 p - 1 ) „ , 2 ,, 

+ ^ , 2 ~ 3 {v + f ) 

l 


L(i + p) (1 + p) 

2 o n ..2 r>2 


/ 


( 1 +p)J p 


Rr- 


3 Rpr)~ O 


A 4 (r j ) = r 5 


2 (1 + p ) 2 P 
3 £R a Fl 2 


e 


(i + p) 

) +r 2 K(i + P )-i] 


Fl 2 R t - 


/ (1 + p) 


2 ’ 


(4.18) 

(4.19) 
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4*(F) = 0, 


(4.20) 


A 2 (T j ) = -r 5 {2 c 2 R a (1 + pf + ft 2 [7 IR£ (1 + p) + 3 R a [p 2 + f) (1 + pf + l 2 Rp (p - 2 )] } 

(4.21) 


x [e+p(/+0] + r4 [ 2i 2^ (1+p) j n 2 [t+p(i+ o] 


2/ 3 (1 + p)' 


2/ 3 (1 + p) 


3 ’ 


(f) 

3 R a pp 2 fl 2 


A 0 (F) = r 5 [e+(/ + e)p] 


whereas the coefficients B n {Tf are defined by 


s 5 (r J ') = r 5 [ 2 c 2 + {rf + e 2 ) ft 2 ] 


ft 2 


+ r 4 


2/3 (1 + p) 
3 R P Q 2 


3 > 


(4.22) 

(4.23) 


2 c 2 ' * V * 

+ T 2 {- (1 + p) 2 [2c 2 + 3ft 2 (p 2 + £ 2 )] + ft 2 [7^ (1 + p) + / 2 (2p - 3)] } 

R/3 


X 


2/3(l + p) z V l 2 
B 4 (r j ) = r 5 


c^y 3^ 2 ft- 

2/3(l + p) 


2 ’ 


3i? Q ft‘ \ 


/ 


+ p 


RpQ 2 

JO- + p) 


, c 2 /? 

B 2 (v j ) = -r 51 u 


/ 3 

ft 2 ^ 
_/(l + p) 


+ ft 2 


i? 3 (F) = 0 , 

7^e , 3R Q ( V 2 +e) 


(4.24) 

(4.25) 

(4.26) 


. 2 / 2 (l + p) 


2/3 


-Rp 3i?pp 


/(1 + p ) 2 2/(1 + p) 2 _ 


Bfr j ) = r 5 ' cRa 


Bo{T j ) = r 5 


/ 2 

3 Rpprfti 2 


(4.27) 

(4.28) 

(4.29) 


_2/ 3 (l + p2)_ • 

The planetoid coordinates are eventually expressed, from the definitions (j4.7[) and (14.81) . in 
the form 

( 1 -P)' 


( 7 2 - r 2 ) + 


(i + p) 


* =l 2 


p = ±/^/7 2 --(7 2 -r 2 + i) 2 . 


(4.30) 


(4.31) 


By numerical analysis of Eqs. fj4.16jl - fj4.17p we have found that, in the Earth-Moon system, 
the only solution where both 7 and T are different from zero is given by 


7 = 0.99999999999996386756, T = 0.99999999999284192083. 


(4.32) 
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These values lead to a tiny departure from the equilateral triangle picture of Newtonian 


theory 

work 


his effect was first predicted in Ref. 


0 ). 


but less pronounced than in our earlier 


151] . where we found a correction of 8.7894 mm on the ^-coordinate and of —4 mm on 


the ^-coordinate [15J. We now find instead, for the planar coordinates of L 4 , 

£(GR) _ £(N) = 2 73 

mm, V (GR) = -1.59 mm. 


(4.33) 


At this stage, we can compare these corrections with those obtained through the method 


outlined by the authors of Ref. 


551 ]. where the position of the Lagrangian points is obtained 


by employing the Einstein-Infeld-Hoffman equation of motion rather than the analysis of 
the zeros of the gradient of the effective potential w in (4.9). As we can see, the correction 
on the ^-coordinate has got the same sign and the same magnitude as the one obtained with 
the method of Ref. 55], while the correction on the 77 -coordinate has got only the same sign, 


because the magnitude is three times bigger. Therefore, it is interesting to note the fact 
that two different methods give exactly the same correction of the ^-coordinate. By taking 
account of Eq. (14.33j) . the resulting values of planetoid distance from Earth and Moon turn 
out to be 

r = 7 / = 3.8439999999998611069 x 10 8 m, s = Tl = 3.8439999999724843437 x 10 8 m. (4.34) 


B. Collinear Lagrangian points 

The position of the collinear Lagrangian points L 4 , L 2 and L 3 is described by the system 
of equations 


dW 
77 = 0, 


= 0 


Following Ref. 15], we know that the vanishing of the / 7 -coordinate implies that 

P 


£ = er — l 


(1 + pY 


( e = ±;L )) 


which in turn leads to the condition 


s = ±(r - el). 


(4.35) 


(4.36) 


(4.37) 


If we substitute relations (I4.36p - (l4.37p into Eq. (I4.12p and initially adopt the choice s = 
(r — el), we obtain an algebraic tenth degree equation where the only unknown is the distance 
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r of the planetoid from the Earth. By setting, as before, r = 7 /, this equation can be written 


down as 


where 


10 


E = 0 , 


c.= 


n =0 

ClO = 1) 

r , _ (V + 4) 

= ~7oTp)’ 

2 c 2 3(7p 2 + 8 p + 2) 


fi 2 / 2 


(1 + pf 


C 7 = 


' e (l _|_ p) {/ 3 Q 2 t 2 ^P + 4) - 3e(l + p)(R a + Rp)] 


a +py 



(4.38) 

(4.39) 

(4.40) 

(4.41) 


(4.42) 


C *6 = e Q 2 p ( 12 (/e — R a — Rp) + p [ 20 /e — 12 (i? Q + Rp)]} 
1 [4p 3 + (l + p)(31p 2 + 14p + l)], 


(1 + P) E 


^ _ ^ c 4 (R q + Rp) + 


/ 5 fi 4 

Rp(18p 2 + 44p + 23) 


(4.43) 


/ 3 fi 2 (l + p ) 2 _ 
3 ep 

(1 + P? 


— 4/e(l + p)(5p + 2) + 3R Q (5p 2 + 12p + 6 ) 

■(7p 2 + 6 p+l), (4.44) 


a = 


2 c 2 

ePQ 2 

1 


[2/(Rp - 2R a ) - e((R a ) 2 + (Rp) 2 )] ^ ] 3 Q 2 ^J + (T+^oj 2 ^ + 2 ^ + ^ - 2e (l + P)\ 


( 1 +py 


:[/e(l + p)(l + 5 p) — 6R q (1 + 2p)] 


P 


(1 + P) : 


:(7p + 3), 


(4.45) 


° 3 = "i«(i + u) + + 6!t)1 “ mi + pfaP 21 ’’ 1 ' 1 + p) 

o 3 

+ 3R a e(5p" — 2p — 1) + Rp[10(l + p) — e(3p" + 44p + 18)]} — ^ ^ 1 (4.46) 

c 2 ( Ac 2 R 1 1 

C ' 2 = J 3 Q 2 | PQ? + 2 ^ + (TTp) 2 ^ 26 ^"^ 2 ~~ + ^p) - (1 + p)]} j ) (4.47) 

2 

Cl = — /30 2 (1 + p) 2 |[ 2c2 -^“(^" + ^ e )(l + P^PQ 2 C 3p“Ra + 2Rp[e(l + p) — (1 + 3p)]|, 

(4.48) 


24 




























Co = 2 


c 2 R a 

IhW 


2 


(4.49) 


whereas in the other case, i.e. s = —(r — el), we end up with the algebraic equation 


10 


J2 d„7” = 0. 


n=0 


with 


D k = C k if k = 10,9,8,0, 

C 2 Rr- 


D 7 = Cj — 6 
Dq = Cq + 24 


03 


W’ 

c 2 i? 


p 


lh = c, - 


2 c 2 R 


m 2 e ’ 

'/3 f 72 o 2 /io „2 , oo. , ml o„2/ 


r [Z 2 fi 2 (18p 2 + 38p + 21) - 2c 2 (1 + p) 2 ], 


d a = c a - 


2c 2 R 


Z 5 h2 4 (l + p) 2 e 
2 c 2 R 


W(l + p) 2 ' 

{4c 2 (l + p ) 2 + 2fl 2 Z 2 [2e(l + p) - ( 6 p 2 + 14p + 9)]}, 


D 3 = C 3 


P 


2 c 2 


+ 


10 


Z 3 ft 2 [1 2 0 2 (1 + p)e (1 + p ) 2 

16c 2 /C 


(3p 2 + 8 p + 6 ) 


d 2 = q 

D V =C\ 


2 2 / 3 12 2 (l + p)’ 

Ac 2 Rpe 


(4.50) 

(4.51) 

(4.52) 

(4.53) 

(4.54) 

(4.55) 

(4.56) 

(4.57) 

(4.58) 


OT(l + p)‘ 

The values Ri {i = 1,2,3) of the distance of the planetoid from the Earth at the libration 
points Li, L 2 , L 3 , respectively, obtained through the solution of Eqs. (14.381) and (14.501) are 
given by 


Rx = 3.2637628817407598555 x 10 8 m, (4.59) 

R 2 = 4.4892056003414800050 x 10 s m, (4.60) 

R 3 = 3.8167471569392170594 x 10 8 m, (4.61) 

whereas the corresponding classical Newtonian values read as 15] 

rx = 3.263762881738878 x 10 8 m, (4.62) 

r 2 = 4.489205600344675 x 10 8 m, (4.63) 

r 3 = 3.816747156939623 x 10 8 m. (4.64) 
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By comparing these values we have 


Rf n) — r[ N 1 = Ri — ri = 0.19 mm at Li, 




(OR) 


R 


(GR) 


r ( N ) _ p 

r 2 — 0,2 


r 2 = —0.32 mm at L 2 , 
= R 3 — r% = —0.04 mm at L 3 . 


(4.65) 

(4.66) 

(4.67) 


Interestingly, the correction on the position of the Lagrangian point L\ is exactly the same 


as the one calculated with the method described in Ref. 


541 ]. where 3 the collinear solutions 


of the three-body problem are studied in the post-Newtonian regime. We believe that, 
according to the definitions involving the ratio of the distances of the planetoid from the 


primaries given in Ref. 


54j , the equations resulting from the application of the method 


developed by the authors of Ref. [54j] (which is the same method used in Ref. 55j) are well 
suited to describe only the position of Li, and the agreement with the corrections presented 
here is a clue supporting our opinion. 


V. QUANTUM EFFECTS ON LAGRANGIAN POINTS 


The analysis of the previous section prepares the ground for a more appropriate definition 
and evaluation of quantum corrections of Lagrangian points, when the underlying classical 


theory of gravity is 
the analysis in Ref. 


instein’s general relativity. For this purpose, we begin by considering 


63] , where the metric tensor components in a co-rotating frame for the 


relativistic restricted planar three-body problem in the post-Newtonian limit were obtained. 
With the notation of our Sec. IV, and coordinates x° = ct , x 1 = £, x 2 = rj, x 3 = (, the result 


3 As shown in Ref. 54[, the general relativity corrections to Li,L 2 may be of order 30 meters in the Sun- 
Jupiter system. However, compared to the Earth-Moon system, a mission to test this effect at Jupiter 
would be exceedingly more expensive and complex to realize and could not even benefit from the use of 
accurate, direct laser ranging from Earth due to the large distance. The effect of the extremely harsh 
Jupiter radiation environment on the test spacecraft (planetoid) should also be considered to evaluate its 
impact on the integrity of the spacecraft and, therefore, the duration of the positioning measurements. 
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in Ref. 


63 ] reads as (cf. Ref. 62]) 


9oo = + + 2 
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2 cg 01 = 1 + 2— + 2-H 2fiij, 

\ r s J 

2 cg 02 = - ( 1 + 2-^ + 2^) 2!2«-8dRy 

#03 = 0 , 


Ra R/3 

p - 


Hi = ~ ( l + 2^ + 2^]iy, i, j = l, 2,3. 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 


The resulting Lagrangian that describes the planetoid motion in the gravitational field of 


Earth and Moon reads as [64J, 


65] 


T 1 ^ da^ dx v 

2^ ^"dTlT 

fi, u=0 


(5.6) 


We now bear in mind that, in light of second line of (4.1), the dimensionless ratio 


u a (r) = — = U m 


(5.7) 


where R a = ^ is the gravitational radius of the primary of mass a, gets replaced by (or 
mapped into) 


V a (r) 


1 + ^2 


1 + ^2 


(h 


(h 


U a (r ) + Ki ^1 + -jj— ) ( U a (r )) 2 + 0(G 3 ) 

^(r) + «i(f/ a (r)) 2 , 


(5.8) 


because the gravitational radius R m of the planetoid or laser ranging test mass is indeed 
much smaller than R n . The same holds for the dimensionless ratio 


£7,(9) = = Uf, 


(5.9) 
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and its effective-gravity counterpart 


Vf(a) 


1 + ^2 


1 + ^2 


{h 


(h 


Up{s) + Ki ^1 + ) (Up(s)Y + 0(G d ), 

Up(s) + Kl (Up(s)) 2 


(5.10) 


By virtue of Eqs. (5.1)-(5.10), we are led to consider the effective-gravity Lagrangian 

Lv = d(l - 2(v„ + V„) - +(f 2 + i; 2 ) + 2 [(VJ 2 + (V s f] 

O i^ a + fy) It 2 i „2 
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{ i + P ) i v ‘ + T?TAT v '- r U R ' v ‘- B ‘ v ' ) 


+ (i + P)~ l7j r 


Rf 3 /-rr \3 Ra \3 


P-nr»(V«r + 


l [ (Ra) 


— -(£ 2 + ?) 2 + C 2 ) 1 + 2 (14 + Vp ) + 1 + 2 ( T 4 + Vg ) 


^ 1 + 2(14 + Vp) 


w 

V 
f l?l 




V(pV a - Vp), 


(5.11) 


(1 +p) 

and the only nontrivial Eulcr-Lagrange equations for the planar restricted three-body prob¬ 
lem are 


Note that, in Refs. 


d / dL v 
dt 1 d£ 


dLv =0 , 4 fife-) -^ = 0. (5.12) 


se 


dt \ di] 


dr] 


13l-ll5|. we have inserted the effective-gravity map (see (5.8) and (5.10)) 


(U a ,Up) ^ (V a ,Vp) 


in the Lagrangian of Newtonian gravity for the restricted planar three-body problem, 
whereas we are here inserting the same map in the Lagrangian of general relativity for the 
restricted three-body problem. The metric tensor with components (5.1)-(5.5) describes, 
within the framework of general relativity, a tiny departure from the Newtonian treatment 
of the restricted planar three-body problem. At that stage, one can recognize that many 
Newtonian-potential terms occur therein; for each of them, we apply the effective-gravity 


map (5.8) and (5.10) to find w' 
Note however that in Ref. 


lat we call a quantum-corrected Lagrangian. 

661 ]. where the authors derive quantum corrections to some 
known exact solutions in general relativity, they find that these metrics differ from the clas¬ 
sical metrics only for an additional term proportional to (Ip) 2 - Within such a framework, 
the running of G at large r has a universal character independent of masses, and there is no 






























room left for K\ in the quantum-corrected Lagrangian. The two schemes are conceptually 
different: quantum corrections to known exact solutions of general relativity do not neces¬ 
sarily have the same nature as quantum corrections of metrics which represent solutions of 
the linearized Einstein equations and which are used in turn to derive equations of motion 
of interacting bodies. The insertion of the map (5.8) and (5.10) in the Lagrangian of general 
relativity for three bodies leads to other terms quadratic in U a and Up, which are of the same 
order of those already present, and hence the resulting Euler-Lagrange equations (5.12) will 
lead to predictions affected by fiq. 

To further clarify this crucial issue we point out that, if we insert the map (5.8) and 
5.10) in the Lagrangian of Newtonian gravity for the restricted planar three-body problem 
rJ, we find, with our notation, the effective potential 


W e g ——(£ 2 + rj 2 ) + c 2 (U a + Up) + Ki((U a y + {Up) 2 ) + 0(G 2 ), (5.13) 

whereas general relativity yields the effective potential (4.9), expressible in the form 
Q 2 r 1 i 

Weff = t (£ 2 + V 2 ) + c 2 [(U a + Up) - -(( U a f + (Up) 2 )\ + 0(G 2 ) 

~ -y(£ 2 + V 2 ) + c 2 (U a + Up) — ~({U a ) 2 + (Up) 2 ) +0(G 2 ), 


(5.14) 


because c u 2 = £(R a + Rp) = 0(G) and, by virtue of (4.4), f2 2 ~ uj 2 + 0(G 2 ). Hence it is 


possib 

Refs. 


13 


o understand why the map (5.8) and (5.10) leads to coordinates of L 4 and L 5 in 


- 15j| pretty close to those of our Sec. IV. 

Moreover, since the work in Refs, lij, 66] has studied three kinds of corrected Newtonian 
potential, i.e. scattering or bound states or one-particle reducible, one has to rewrite the 
map (5.8) and (5.10) in the form 


V a ~u a + U + 7 -K^j ( U a ) 2 + 0(G 2 ), 
Vp ~ Up + + 2^°^ (^(s) 2 + 0(*^ 2 )> 


(5.15) 

(5.16) 


which, upon defining K\ = k\ + -k 0 , takes the form 

V q ~G q + ki(G q ) 2 + 0(G 2 ), 


(5.17) 


Vp-Up + K^Upf + OiG 2 ), 


(5.18) 
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where 


k 2 


(lp) 2 TT _Cl(r 2 \ r ^ P ) 

- 5 — Ua — 0(G J,/v2 - 7 


-Up = 0(G 2 ), 


and the parameter Ko vanishes in the scattering and one-particle reducible cases 66 ], whereas 
it equals —1 for bound states [12]. Remarkably, since =3, and k 0 = — 1 for bound states 
[eJ ]. this simple calculation shows that the insertion of the map (5.15) and (5.16) into the 
Lagrangian of Newtonian gravity for the three-body problem leads to the effective potential 


Weff = y (£ 2 + V 2 ) + c 2 [{U a + Up) - ^(( U a ) 2 + (t^) 2 )] + 0(G 2 ), (5.19) 


which has the first three terms in common with the effective potential (5.14) of general 
relativity. As far as we can see, this is evidence in favour of inserting the effective-gravity 
map into the Lagrangian, and in favour of considering the values of K\ and k 0 appropriate 
for bound states both below and in the analysis of Secs. II and III. Now we set to zero all 
time derivatives of £ and 77 in Eqs. (5.12), we define the real numbers 7 and T as in Sec. IV 
and solve numerically the resulting algebraic system for such numbers. This method leads 
to the following values: 


• noncollinear Lagrangian points: 

The planar coordinates of equilibrium points L 4 and L 5 read as 

£4 = £5 = 1.8752814880352634039 x 10 8 m, (5.20) 

774,5 = ±3.3290016521227759284 x 10 8 m, (5.21) 

which means that the differences with respect to the corresponding values provided 
by the Einstein theory, that plays in this scheme the role of the classical theory of 
reference, read as (cf. Eq. (4.33)) 

£4 — = —1.46 mm £4 — £^ A 1 = 1.27 mm, (5.22) 

774 — 774 GR) = —0.86 mm 774 — 77 ^ = —2.45 mm. (5.23) 


• collinear Lagrangian points : 

For the libration points L 1 , L 2 and L 3 , respectively, we have found that 

R\ = 3.2637628817345938976 x 10 8 m, (5.24) 
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R' 2 = 4.4892056003375634274 x 10 8 m, (5.25) 

R! 3 = 3.81674715692440418189 x 10 s m. (5.26) 

Thus, bearing in mind Eqs. (I4.59p ~ fj4.61j) . the differences with respect to the values 
expected from general relativity are 


R[ — R,\ — —0.62 mm, 
R' 2 — R 2 = —0.39 mm, 


R 3 — R 3 = —1.48 mm. 


(5.27) 

(5.28) 

(5.29) 


Another important issue concerning both noncollinear and collinear Lagrangian points, 
consists in the fact that we have checked numerically that the corrections (I5.22p . (15.231) and 
(I5.27p - (l5.29p do not change if we set k 2 = 0 in the Euler-Lagrange equations (15. 12j) . because 
k 2 weighs the dimensionless ratios and which are extremely small at large values 
of r and s. 

VI. GENERAL RELATIVITY VS. EFFECTIVE GRAVITY: CONCLUDING RE¬ 
MARKS AND OPEN PROBLEMS 


The fir st p art of this paper contributes to make more realistic the model outlined in 


Refs. 


13H15| by considering the gravitational presence of the Sun as a perturbing effect 


for the Earth-Moon system. In fact we have shown that also in the quantum regime the 
presence of the Sun makes the planetoid ultimately escape from the triangular libration 
points, which therefore can be considered as “stable” equilibrium points only during the 
length of observations. Unless we consider solar radiation pressure, from Eqs. (12.181) (I2.20j) 
we have obtained a plot describing the spacecraft motion about L 4 (Fig. [3]) which is slightly 
modified if compared with the corresponding classical one (Fig. [[]). If we instead take into 
account the solar radiation pressure (13.11) . the differences between classical and quantum 
theory become more evident. The presence of solar radiation pressure in the classical case, 
in fact, makes just the planetoid go away from the Lagrangian points L 4 more rapidly (see 
Fig. [5]), but in the quantum case, before escaping away from the libration point L 4 , the 
planetoid is characterized by a less chaotic and irregular motion, as is clear from Fig. [61 
This feature remains true also if we consider several initial velocities for the planetoid (Figs. 
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[7] and IHJ) . In particular, we have shown that the reduction of the envelope of the planetoid 
motion becomes more evident in the quantum case. After that, we have calculated the 
impulse needed for the stability of the spacecraft at L 4 both in the classical and in the 
quantum regime. These two values, as witnessed by Eq. (I3.10p . are a little bit different and 
therefore they suggest sending two satellites at L 4 and L 5 , respectively, and checking which 
is the impulse truly needed for stability, in order to find out which is, between the classical 
and the quantum one, the best theory suited to describe these phenomena. 

In the second part of the paper, we first perform a comparison between Newtonian gravity 
and general relativity, since of course the latter is the most successful theory describing 
gravitational interactions, at least in the solar system. By evaluating the points where the 
gradient of the potential (14.91) vanishes, we have solved the algebraic equation describing 
the position of Lagrangian points. The distances of noncollinear Lagrangian points from the 
primaries are given in terms of the solutions of Eqs. (14.14p (j4.15[) (or equivalently (I4.16p - 
(I4.17P ). As is clear from Eq. (14.331) . we have obtained corrections of the planar coordinates 
of the triangular libration points of the order of few millimeters and, in particular, the 
correction on the ^-coordinate is exactly the same as the one obtained through the method 


outlined by the authors of Ref. 


551 ]. As far as collinear Lagrangian points are concerned, 


we have to focus on Eqs. (14.381) and (I4.50p . from which we have evaluated corrections of 
the distances of the planetoid from the Earth of the order of few meters. In particular, the 
correction concerning L\ is exactly the same as the one obtainable with the method of Ref. 

54]. 

In the last part of our paper, we have outlined the features of a quantum theory whose 


underlying classical theory is represented by general relativity and not, as before, by New¬ 
tonian gravity 


I 21 15]. I 11 other words, we have dealt with a theory involving quantum 
corrections to Einstein gravity, rather than to Newtonian gravity. In fact, by applying the 
map (15.81) and (15.101) to the Lagrangian 


that general relativity provides for the re¬ 
stricted three-body problem, we have ended up with the quantum corrected Lagrangian 
(15. lip which, by the means of Euler-Lagrange equations (15.121) together with the conditions 

£ = ^ = C = £ = ? 7 = C = C = 0) has led us to the few millimeters corrections (I5.22p . (I5.23P 

41 

and (I5.27p - (l5.29p . obtained by using the quantum coefficients fiq = —1/2 and k 2 = tt— of 


the bound-states potential discussed in Ref. 


107T 


12], We stress that, within this new scheme, 


quantum corrections on Newtonian quantities have been obtained through the algebraic sum 
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of quantum corrections to general relativity (obtained by (15.121) 1 and general relativity cor¬ 
rections to Newton’s theory (obtained by the solutions of Eqs. (I4.16p . (I4.17[) . (I4.38jl and 


(fmi l. In ot 


in Refs. 


13 


rer words, in this new approach we are no longer using the method developed 


15 ]. where the map (5.8) and (5.10) was only inserted into the Newtonian 


Lagrangian of the restricted three-body problem, for simplicity. The possibility of mapping 
the effective potential of Newtonian gravity into an effective potential similar to the one of 
general relativity (cf. (5.19) and (5.14)) adds evidence in favor of the choice of and k 2 
appropriate for bound states 12], We also believe it is important to stress that we have 
used modern packages for dealing with all the coupled algebraic equations presented in this 
paper, verifying eventually that the putative solution does satisfy the original set. 

In conclusion, as far as we can see, the implications of calculations presented here and in 
our previous work are as follows. 


A. Noncollinear Lagrangian points 

(i) Quantum corrections to Newtonian planar coordinates of L 4 and L 5 and to the cor¬ 
responding general relativity values are (1.27 mm, -2.45 mm) and (-1.46 mm, 0.86 mm), 
respectively. If it were possible to obtain an experimental verification of this prediction, 
it might provide encouraging (but not conclusive) evidence in favor of effective theories of 
gravity with values Ki = —1/2 and k 2 = 41 /(107r) of the bound-states potential. 

(ii) If the departure from the equilateral triangle picture of Newtonian theory is instead close 
to (2.73 mm, -1.59 mm) for the planar coordinates of L 4 and L 5 , this may provide another 
nontrivial test of classical general relativity. 


B. Collinear Lagrangian points 

(iii) If the differences with respect to the classical Newtonian distances of the planetoid at 
Li, L 2 and L 3 from the Earth were of order (0.19 mm, -0.32 mm, -0.04 mm), we would have 
to consider this fact as another confirmation of general relativity. 

(iv) In the case in which the quantum theory is ruled by the Lagrangian (15.lip , the quantum 
corrections to general relativity become (-0.62 mm, -0.39 mm, -1.48 mm) at Li,L 2l L 3 , 
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TABLE I: Distances r 4 from the Earth and planar coordinates of the planetoid at 

all Lagrangian points Lj in the classical Newtonian theory, General Relativity and the 
quantum regime, the latter being obtained through the Lagrangian Ly (15.111) and with the 
coefficients of the bound-states potential kj = —1/2 and H 2 = 41 /(107r) taken from Ref. 

Ir¬ 


Li 

Newtonian Gravity 

General Relativity 

Quantum regime 


r\ = 3.263762881738878 x 10 s m 

n = 3.263762881740760 x 10 s m 

ri = 3.263762881734594 x 10 8 m 

Li 

£1 = 3.217044369761366 x 10 s m 

£1 = 3.217044369763247 x 10 8 m 

£1 = 3.217044369757081 x 10 8 m 


II 

o 

O 

II 

£ 

O 

II 

s? 


r 2 = 4.489205600344675 x 10 8 m 

r 2 = 4.489205600341480 x 10 8 m 

r 2 = 4.489205600337563 x 10 8 m 

l 2 

6 = 4.442487088367163 x 10 s m 

£2 = 4.442487088363968 x 10 8 m 

£ 2 = 4.442487088360051 x 10 8 m 


V2 = 0 

V2 = 0 

m = 0 


r 3 = 3.816747156939623 x 10 s m 

r 3 = 3.816747156939217 x 10 8 m 

r 3 = 3.816747156924404 x 10 8 m 

l 3 

$3 = -3.863465668917136 x 10 s m 

£3 = -3.863465668916729 x 10 8 m 

£3 = -3.863465668901917 x 10 s m 


??3 = 0 

T)3 =0 

V 3 = 0 


r 4 = 3.844000000000000 x 10 8 m 

r 4 = 3.843999999999861 x 10 8 m 

ta = 3.843999999985078 x 10 8 m 

La 

£4 = 1.875281488022488 x 10 s m 

£ 4 = 1.875281488049864 x 10 8 m 

£ 4 = 1.875281488035263 x 10 8 m 


r/A = 3.329001652147382 x 10 8 m 

r/A = 3.329001652131416 x 10 8 m 

r/A = 3.329001652122776 x 10 8 m 


r 5 = 3.844000000000000 x 10 8 m 

r 5 = 3.843999999999861 x 10 8 m 

r 5 = 3.843999999985078 x 10 8 m 


$5 = 1.875281488022488 x 10 s m 

£ 5 = 1.875281488049864 x 10 8 m 

£5 = 1.875281488035263 x 10 8 m 


V5 = -3.329001652147382 x 10 s m 

7/5 = -3.329001652131416 x 10 s m 

775 = -3.329001652122776 x 10 s m 


respectively. 

A summary of all quantities involved in this paper along with all corrections discussed 
above is reported for clarity in Tabs. [T] and [Til As far as we can see, our detailed calculations 
show clearly that the measurement we are proposing in the Earth-Moon system represents 
a new testbed for general relativity and effective field theories of gravity. 
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TABLE II: Corrections to the distances of L i, L 2 and L 3 from the Earth and to the planar 
coordinates of L 4 and L 5 . In the first column we have the general relativity corrections to 
Newtonian theory as obtained from Eqs. (I4.38E (I4.50L (I4.16ji and (I4.17L respectively. The 
second column shows quantum corrections to general relativity given by the 
Euler-Lagrange equations involving the Lagrangian Ly (15. lip and the quantum coefficients 
of the bound-states potential. The last column displays quantum corrections to Newtonian 
values calculated as the algebraic sum of the corresponding quantities in the two previous 

columns. 


U 

General Relativity-Newton 

Quantum-General Relativity 

Quantum-Newton 

Li 

0.19 mm 

-0.62 mm 

-0.43 mm 

L2 

-0.32 mm 

-0.39 mm 

-0.71 mm 

L3 

-0.04 mm 

-1.48 mm 

-1.52 mm 

l 4 

(2.73 mm, -1.59 mm) 

(-1.46 mm, -0.86mm) 

(1.27 mm, -2.45 mm) 

L5 

(2.73 mm, -1.59 mm) 

(-1.46 mm, -0.86mm) 

(1.27 mm, -2.45 mm) 


Appendix A: Evaluation of the gradient of W in general relativity 


By virtue of formulas (4.7) and (4.8), one finds 

d 


pi 


- {r - n = -rr-v-1 [i+ {1 + p)J , 


d_ 

dp 

d_ 

dp 


V p ) = — pr p 2 , 


s P) = -p S p 2 , 


(Al) 

(A 2 ) 

(A3) 

(A4) 


that we have computed with p — 1, 2, 3 to obtain Eqs. (4.12) and (4.13) for the components 
of grad(hL), where the functions W1...W4 are defined by 

^ V 


W^,p,r) = + 


Q 4 ^ 2 + p 2 ) ^R n 7 R«m 2 1 


2 c 2 


+ 3fi 2 e— + n 


r 2 (1 + p) r 




pi 


— / 2 /? f — 

(1 + p) J r 3 { a \ r 


l)+fl 2 -^R a (e + P 2 ) 


7 Rpl£ + 3 


p Rpl 2 p 2 , (2 - p) 2 

-Kpl 


+ 


2 (1 + p) 2 (l + p ) 2 r 2 2(1 + p ) 2 


(AS) 
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W 2 (£,ri,s) = 


Rn 7 RJQ 2 1 


3Q 2 ^ - ' ** 


e- 


s 2 (1 + p) s 

(T h) (t 

2„2 


7 TfyZf 3 1 RpFri- 

2 (1 + p) + 2 (1 + p) 2 s 2 


- lj + ft 2 

, (2p-l) 

2(1 + p) 2 


-^Me+v 2 


Rsl 2 


(A6) 


W 3 {£,rj,r) = 


r ^2 ^ / a2 9\ ^ -^q; -^q 

n + 5? ( « +J ' ) + — (t 

3q2^2 _|_ _ "fR.rjlfiT 1 


- 1 


30' 


77 


2 (1 + p) r 3 


77/jZ 2 ^ 2 p f n 

2(1 + p) 2 r 3 \ r 


2 + 


77^ 2 fi 2 (2 — p) 
2(1 + p) 2 r 3 


(A7) 




c 2 Rb 


Rpi 2 n 2 i 


^ - 1) + - jW + rj*)- 

s J s 2 s 3 

_ 2 V R ^ 2 (2p “ 1} 


7?^ 7 7y£0 2 1 


2 (1 + p) s 3 


2(1 + p) 2 s 3 V « 


2(1+ p) 2 s 3 
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